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Abstract. Transport properties of two coupled Josephson junctions driven by ac 
currents and thermal fluctuations are studied with the purpose of determining dc 
voltage characteristics. It is a physical realization of directed transport induced by 

O a non-biased zero averaged external signal. The ac current is applied either to (A) 

only one junction as a biharmonic current or (B) is split into two simple harmonic 

! K components and separately applied to respective junctions. We identify regimes where 

junctions can operate with the same as well as opposite signs of voltages. A general 

^ observation is that in the same parameters regimes, the scenario (B) is more efficient 

"£i in the sense that the induced dc voltages take greater values. 
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1. Introduction 

Noisy transport in periodic arragements pQ is widely present in many nowadays branches 
of science - in physics, biology, chemistry, economy and many others. On the physical 
ground, the periodicity itself can be associated either with space degrees of freedom like 
in the crystals, optical lattices, systems of ring topologies or with time-periodic drivings 
like ac currents, magnetic or electric fields, rocking and pulsating forces to name but 
the few. It also can be present in these both domains. Typical realizations can range 
from biophysics [2] with the description of biomotors movement on asymmetric periodic 
microtubules [3] or transport inside ion channels [4], to the present experiments with 
optical lattices |5j|6], quantum mesorings [7] or Josepshon junctions [8]. 

The Josephson effect is known for a half of the century [9]. Since that time it 
has been utilized for the definition of the voltage standard |10j or for more practical 
devices as elements in high speed circuits [8] or even for the future applications in 
quantum computing devices [Uj. Surprisingly, after 50 years of intensive theoretical and 
experimental research, we are still able to find new and uncommon phenomena even in a 
simple system of two weakly connected superconductors. Recently, the counterintuitive 
phenomenon of absolute negative conductance (ANC) has been reported in the single 
driven, resistively, and capacitively shunted Josephson junction device subjected to both 
a time-periodic (ac) and a constant biasing (dc) current [12]. The ANC phenomenon 
has been confirmed by the suitable experiment with a Josephson junction setup [13] and 
very recently with ultracold atoms in optical lattices [5J. Other aspects of anomalous 
transport phenomena like the occurrence of a negative differential conductance and the 
emergence of a negative nonlinear conductance in the nonequilibrium response regime 
remote from zero dc bias have been studied in a series of papers jH] . The influence of the 
unbiased biharmonic ac current on a single junction has been considered in Refs. [15J. 
Recently the dynamics of the phase difference of coupled junctions has been addressed 

This work is organized as follows. In section [2j we present the model of two 
interacting junctions. Next, in section [3j the numerical investigation of transport 
properties for two scenarios (A) and (B) of drivings applied to two coupled junctions is 
compared. The paper ends with summary and conclusions in section 4. 

2. Model of driven interacting junctions 

From a more general point of wiev, we explore the system consisting of two subunits 
(subsystems) interacting with each other. The system is driven out of its equilibrium 
state by an external force. As a particular realization of this idea we propose two 
resistively shunted Josephson junction devices characterized by the critical Josephson 
supercurrents (I cl ,I c2 ), resistances (i?i, R 2 ) and phases (0i, <p 2 ) [IB]- A schematic circuit 
representing the model is shown in figure [T] The system is externally shunted by the 
resistance R 3 and driven by two current sources I\(t) and I 2 {t) acting on the first and 
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second junction, respectively. We also include into the model Johnson-Nyquist thermal 
noise sources £i(t),£, 2 (t) and ^(t) associated with the corresponding resistances Ri,R 2 
and i?3 according to the fluctuation-dissipation theorem. 
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Figure 1. The system of two coupled Josephson junctions characterized by the 
critical Josephson supercurrents (J c i, I C 2), resistances (R±, R2), shunted by the external 
resistance R3, influenced by Johnson-Nyquist thermal noise sources (^,i{t) , ^(t) , ^(t)) 
and driven by the external currents (Ii(t), ^(t)). 



The beauty of the standard Josephson theory lies in the simplicity of the model. In 
the semiclassical regime, when the spatial dependence of characteristics can be neglected 
and photon-assisted tunnelling phenomena do not contribute, the so-called Stewart- 
McCumber model [19] holds true (for the extensive discussion on the validity of the 
model we refer the reader to |10|). In this regime, one can use the classical Kirchhoff 
current and voltage laws, and two Josephson relations to derive two evolution equations 
for the phases 0i = 0i(£) and 02 = 02 (t). The dimensional version of the equations is 
presented in Ref. {17] . Here, we recall their dimensionless form, namely, 

0i = 2i(r) -/ c isin0i + a[I 2 (r) -/ c2 sin0 2 ] + VDr]i(r), (la) 

02 = afi[I 2 {r) - I C 2 sin 2 ] + a[ii (r ) - I cl sin 0i] + ^a/3D r] 2 (r) , (16) 

where 0j = 0j(t) for i = 1,2 and the dot denotes a derivative with the respect to the 
dimensionless time r expressed by the dimensional time t as 
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All dimensionless currents Ii(t),I 2 (t),I c i and I c2 are in units of I c , e.g. I c \ — > I c \jl c . 
Thermal equilibrium noise sources related to the resistances Ri,R 2 , R3 are modelled here 
by the independent ^-correlated zero-mean Gaussian white noises £j(£) (i = 1, 2, 3), i.e., 
(t,i{t)£,ji s )) = $ijfi{t ~ s ) f° r hj ^ {1; 2, 3}. The straightforward assumption of identical 
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temperature T felt by all parts of the set-up allows for the reduction of the number 
of original noises ^1,^2^3 (see figure II]), to their linear combination 771 and r/ 2 in the 



equations (la) and (lb). 



The reader would find easier to understand this scenario within a pure mechanical 
picture. The dynamics of the phase difference can be mapped onto the motion of the 
Brownian particle. In this mechanical analog the correspondence between position X\ 
of the first particle with the phase difference fa of the first junction can be settled and 
the position x 2 of the second particle can mimic the phase difference 02 of the second 
junction. If we imagine two interacting particles moving along the periodic structure 
than the most significant quantifiers describing their transport properties would be the 
average velocities of the first v\ = (fa) and second v% = (02 ) particle, respectively. 
In terms of the Josephson junction system it corresponds to the dimensionless long- 
time averaged voltages V\ = (fa) and v 2 = (fa) across the first and second junctions, 
respectively (from the Josephson relation, the dimensional voltage V = (h/2e)d<j)/dt and 
therefore d<j)/dr = V/Vq). The junction resistances (or conductance) translates then into 
the particles mobility. Moreover, the phase space of the deterministic system is three- 
dimensional {0i (r), 02 (t), ut} and therefore give rise to possible chaotic evolution which 
is the key feature for anomalous transport [T2| HU [20] . 

2. 1 . Identical junctions 

Without loss of generality, we can reduce a number of parameters assuming that two 
junctions are identical with R\ = R 2 and I c \ = I c2 = 1. In such a case a(3 = 1 and 



equations (la) and (lb) take symmetric form 

0i = Ii(t) - sin 0i + a[I 2 (r) - sin0 2 ] + VI) r?i(r), (4a) 

02 = h(r) -sin0 2 + a[Ji(r) -sin0i] + VDt] 2 (t). (Ab) 

The parameter a = #2/(^2+^3) £ [0, 1] plays the role of coupling strength between the 
junctions and can be tuned by the variation of the external resistance R 3 . When a = 
the set of equations (4) decouple into two independent equations. It can be realized 
taking R 3 — > 00. The opposite situation with two fully coupled junctions can be worked 
out by designating R3 = 0. The noise strength D can be tuned by temperature. The 
currents Ii(r) and -^(t) are energy sources pumped into the system and can be applied 
to one or to both junctions. 

2.2. External current driving 

The trivial way to induce the dc voltage across both junctions is to apply the dc current 
to both junctions separately (in the mechanical analog, it corresponds to the static 
force). It seems that we can also do it by applying the dc current to one junction only 
but have to make sure that coupling is strong enough to call out the response on the 
other junction too. This, however, seems to be rather uninteresting and a well known 
solution. What if we abandon simple intuitive possibilities? We can exploit the well 
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known ratchet effect [2T] and induce the non-zero dc voltage by applying a zero-mean 
external current. We consider two scenarions. In the first scenario (A), the ac driving 
to applied to only one of the junctions [151 [T7] . namely, 

hij) = ai cos(wr) + a 2 cos(kuT + 9), I 2 (t) = 0. (5) 

where 9 is the relative phase between the driving currents and A; is a real number. 

In the second scenario (B), the external current is split into two simple harmonic 
components applied to two respective junctions, namely, 

hij) = ai cos(wr), hij) = a 2 cos(kcoT + 9). (6) 

We know that the symmetric driving cannot itself induce the non-zero dc voltage. 
However, we expect that the coupling between junctions would have to play the crucial 
role in the dynamics of the total system and a non-zero dc voltage could be generated 
for a > 0. We ask which of two scenarios (5) or (6) is more efficient in the sense that the 
induced dc voltages have greater amplitudes. In the method (5) we have the possibility 
to induce the non-zero dc voltage just by the ratchet effect, cf. the detailed discussion 
in Ref. [17J. In this case, even for a = we still can find non-zero dc voltage across 
the first junction. In the scenario (6) the separated symmetric ac currents cannot alone 
induce non-zero voltage in the decoupled junctions. Setting the parameter a ^ we 
effectively incorporate the ratchet effect and in turn create the prospect of dc transport 
in the system. 

3. Dc voltage characteristics 

Stochastic differential equations (4) cannot be handled by known analytical methods. 
For this reason we have carried out extensive numerical simulations. We have used the 
2 nd order Stochastic Runge-Kutta algorithm with the time step of about 10~ 3 • {2tt/uj). 
The initial phases 0i(O) and 02(0) have been randomly chosen from the interval [0, 271"]. 
Averaging was performed over 10 3 — 10 6 different realizations and over one period of 
the external driving 2ir/u. Numerical simulations have been carried out using CUDA 
environment on desktop computing processor NVIDIA GeForce GTX 285. This gave 
us possibility to speed up the numerical calculations up to few hundreds times more 
than on typical modern CPUs [22]. Below, we present results for a fixed frequency 
multiplier k = 2. This will reflect the typical biharmonic driving studied previously 
for Hamiltonian systems [23], systems in the overdamped regime [24"1 123] and for the 
moderate damping [261 127] . If the given parameter is not addressed directly in the 
plot we will keep the constant values as follows: the noise strength (or equivalently the 
dimensionless temperature) D = 0.001, the frequency of the ac driving uj = 0.03944, the 
coupling strength a = 0.56, the relative phase 9 = tt/2 and the amplitudes a\ = a 2 = 1. 
In the long time limit, the averaged voltages (</>j(r)) can be presented in the form 
of a series of all possible harmonics, namely, 

oo 

Jim (4>i{r)) = Vi + Y^ Vi(nurr), i = l,2, (7) 

n=\ 
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Figure 2. (color online) The stationary averaged dc voltages v\ and V2 across the 
first and second junction. The dependence on the external ac current amplitudes 
ai and a-i are presented in panels (a) and (b) for the driving (5) acting on the first 
junction only and in panels (c) and (d) for the driving (6) split between two junctions. 
Other parameters read: the dimcnsionless temperature D — 0.001, the frequency 
ui = 0.03944, coupling strength a = 0.56, the relative phase 9 = tt/2 and the frequency 
multiplier k = 2. 



where V{ is a dc (time- independent) component and Vi(nut) are time-periodic functions 
of zero average over a basic period. For high frequency u (i.e. fast alternating currents) 
the averaged dc voltages are zero: very fast positive and negative changes of the 
driving current cannot induce the dc voltage and only multi-harmonic components of 
the voltages can survive. In addition, if both amplitudes ai,a 2 are smaller than the 
critical supercurrents, from the structure of the model (4) it follows that the net voltage 
will be zero or very close to zero. 

In figure [2] the long-time averaged dc voltages across the first (t>i) and second (^2) 
junctions are shown in the amplitudes parameter plane {01,02}. There is clearly zero 
average voltage for small values of both amplitudes. However, for larger amplitudes in 
both scenarios (5) and (6) we can recognise four operating regimes where: 

(i) v 1 > and t> 2 > 0, 
(ii) v 1 < and t> 2 < 0, 

V\ < and t> 2 > 0, 

Vi > and t> 2 < 0, 
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Figure 3. (color online) The stationary averaged dc voltages v\ and i>2 across the first 

( blue line) and second ( red line) junction, respectively. The dependence 

on the external ac current amplitudes a,\ and a 2 are presented for the driving (5) in 
panels (a) and (b) and for the driving (6) in panels (c) and (d). Other parameters 
read: the dimensionless temperature D = 0.001, the frequency to = 0.03944, coupling 
strength a — 0.56, the relative phase = ir/2 and the frequency multiplier k — 2 and 
the amplitudes e^ — a 2 — 1. 



Of course, the quantitative picture is different. One can easily see that for the case 
(5), transport properties of the first (driven) junction have more complicated strips-like 
structure with larger area of negative voltage. For the second (non-driven) junction, 
the dc voltage v<i can be twice - three times greater than V\. On the other hand, in 
the case (6), the regimes of negative voltage are smaller. We emphasize that such 
complicated regimes of islands and tongues of negative and positive dc voltages are not 
just rare occurrences: they can be verified with numerically arbitrarily-high-accuracy 
calculations and over extended intervals in the parameter space. 

In figure |3l we have selected the region where for the scenario (5) the dc voltage on 
the first (driven) junction stays negative or zero (more precisely, so small to be negligible) 
throughout the presented range of both amplitudes. The second (non-driven) junction 
shows all possible working states with the negative, positive and zero dc voltage (see 
left panels of figure [3]). In the same region but for the second scenario (6), the first 
junction driven by the current cosur stays positive for all presented values of both 
current amplitudes Oi and 02, while the junction driven by the current cos(2o;r + 7r/2) 
can assume positive and negative values. What is also striking is that in the case (6) 
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Figure 4. (color online) The stationary averaged dc voltages v\ and t>2 across the first 

( blue line) and second ( red line) junction, respectively. The dependence on 

the noise strength (or dimensionless temperature) D (upper panels) and the coupling 
strength a (bottom panels) are presented for the driving (5) in panels (a) and (b), 
and for the driving (6) in panels (c) and (d). Other parameters if not addressed 
directly in the plots read: the dimensionless temperature D = 0.001, the frequency 
oo = 0.03944, coupling strength a — 0.56, the relative phase 9 = tt/2, the current 
amplitudes a,\ = a% = 1 and the frequency multiplier k = 2. 



the voltage characteristics of both junctions change with some synchrony: the voltages 
increase or decrease when one of the current amplitude varies. This feature, in turn, is 
not found in the scheme (5). 

As the next point of analysis, we ask about the role of thermal fluctuations. It is 
presented in upper panels of figure |4j In this regime we can note the voltage reversal 
across the second junction: the voltage v 2 can change its sign from negative to positive 
values when temperature is increased. On the other hand, the voltage v\ is always 
negative for the scenario (5) and is always positive for the scenario (6). For high 
temperature, both voltages tend to zero. Next, we address the issue of whether, and 
to which extent, the coupling strength a can influence voltage properties. The results 
are depicted in bottom panels of figure |4j The first note is non-monotonic and irregular 
dependence of both volatges on a with several minima and maxima. In the scenario 
(6), a step-like dependence of the voltage across the first junction is observed for small 
values of the coupling a. 
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4. Summary 

With this study we numerically analyzed the role of ac current drivings on transport 
properties of the resistively shunted two coupled Josephson junctions. We identified a 
rich variety of the dc voltage characteristics in the parameter space where transport can 
be experimentally monitored. We have detected regions displaying positive and negative 
dc voltages, which form complicated structures in the parameter space. We have mainly 
concentrated the analysis on impact of selected regimes in the parameter space on voltage 
properties. Other regimes of parameters also modify voltage characteristics but here we 
do not present all varieties. A general observation is that in the same parameters 
regimes, the biharmonic ac driving applied only to one junction results is a smaller dc 
voltage than in the case when the ac current is split into two simple harmonics and each 
applied to respective junctions. 
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